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Abstract Open-access microcavities are emerging as a new 
approach to confine and engineer light at mode volumes down 
to the A 3 regime. They offer direct access to a highly confined 
electromagnetic field while maintaining tunability of the system 
and flexibility for coupling to a range of matter systems. This 
article presents a study of coupled cavities, for which the sub¬ 
strates are produced using Focused Ion Beam milling. Based 
on experimental and theoretical investigation the engineering of 
the coupling between two microcavities with radius of curvature 
of 6pm is demonstrated. Details are provided by studying the 
evolution of spectral, spatial and polarisation properties through 
the transition from isolated to coupled cavities. Normal mode 
splittings up to 20 meV are observed for total mode volumes 
around 10 x A 3 . This work is of importance for future devel¬ 
opment of lab-on-a-chip sensors and photonic open-access 
devices ranging from polariton systems to quantum simulators. 



Spectral engineering of coupled open-access microcavities 

L.C. Flatten* *, A.A.P. Trichet, J.M. Smith 


1. Introduction 

Optical microcavities, owing to their small confinement vol¬ 
ume, serve to enhance the light-matter interaction. They 
significantly improve the control and sensing capacities 
for electronic systems situated in the confined region. Epi¬ 
taxial and monolithic microcavities have made substantial 
contributions to various fields such as quantum optics [1], 
lasing [2], optomechanics [3] and optical sensing [4]. Build¬ 
ing on improvements in the fabrication methods of these 
devices, it has been possible to move beyond the single 
cavity design towards coupled cavities [5,6] and on-chip 
photonic circuits [7-9]. As striking examples, coupled cav¬ 
ities have been used to realise highly efficient sources of 
entangled photon pairs based on semiconductor quantum 
dots [10], macroscopic Josephson oscillations of exciton 
polaritons [11] and photonic graphene [12,13]. They are 
a versatile building block for the realisation of enhanced 
sensors [14], single photon sources [15], effective hamilto- 
nian engineering [16] and may allow advances in photonic 
quantum simulations [17], the study of Tonks-Girardeau 
gases [18] and many-body physics in general [19]. 

Much of the focus for microcavity research has been 
on monolithic solid state structures, such as micropillars 
grown by epitaxial methods. Recently however, an interest 
has emerged in open-access microcavities because of their 
facile tunability and flexibility to couple a range of elec¬ 
tronic systems to the photonic mode [20-27]. Open cavities 
bring the additional advantage over micropillar structures 
that confinement of the optical mode can be achieved on 


sub-micron scales without either disrupting the integrity 
of the electronic system or causing scattering losses that 
inevitably occur with conventional etching and patterning 
techniques. Adding to that, a high degree of control over 
the topography [28], low mode volumes down to A 3 and 
high finesses up to 2 x 10 5 were achieved with the current 
fabrication techniques [29-31]. Recent demonstrations of 
polariton formation with quantum wells [32] and Purcell 
enhancements with 2D materials [33] have evidenced the 
maturation of the open-access microcavity design. 

In this work, we describe the first study of the engi¬ 
neered coupling between two of such open cavities forming 
a photonic molecule [5,6,8]. We present the experimental 
formation of strongly coupled ‘supermodes’ in the open- 
access geometry and model the results using a modemixing 
formalism. Two identical cavities are gradually merged mov¬ 
ing from two uncoupled oscillators (features fully separated) 
to a single one (features fully merged). In between, the hy¬ 
bridisation of the fundamental mode as well as the first 
four transverse excited states is analysed. Moreover, taking 
advantage of the intrinsic tunability, we demonstrate the 
anti-crossing behaviour of this strongly coupled system by 
simply changing the angle between the two mirrors. This 
study paves the way to explore effects such as the uncon¬ 
ventional polariton blockade [34-36] and create advanced 
structures such as rings and coupled arrays for future devel¬ 
opments in the field of open-access nanophotonics. 
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Figure 1: a) FIB micrograph of coupled cavities templates. 
From left to the right the distance d between the two features 
is changed from d = (10,0) pm over (0,0) pm to (0,10) pm. 
b) Experimental setup showing control over cavity length 
L and angle 6 between the two mirrors, c) Length initiali¬ 
sation routine with white light transmission, d) Laser trans¬ 
mission and equivalence between length and spectral sweep. 


2. Chip design and experimental setup 


shape with nanometric accuracy [28]. These templates are 
then coated with Distributed Bragg Reflectors (DBR’s) at 
the Thin Film Facility of the Department of Physics (Univer¬ 
sity of Oxford) via sputter deposition. In this work, we used 
10 pairs of Si 02 /Ti 02 giving a macro-reflectivity above 
99.7% at the center wavelength of A = 637 nm, correspond¬ 
ing to a maximum achievable finesse ^ 1000. 

In order to study the gradual coupling behavior between 
two cavities, we produced chips with 13x13 features with 
varying inter-cavity distance d. The single feature shape 
was chosen to be the isophase surface of a Gaussian TEMoo 
mode function (for further information on the exact shape 
and cavity merging parameters see Supplement 1). Lig. SI 
a) displays a FIB micrograph of a selection of these features 
on a silicon substrate, silica being non-conductive. In this 
example, the parameter d is scanned from d = 10 pm along 
the x axis (left cavity) to d = 10 pm along the y axis (right 
cavity), going through fully merged cavity for d = 0 pm 
along both axes. The spectral properties of the fully merged 
cavities and separated ones are expected to be similar since 
they feature the same depth (400 nm) and radius of curvature 
R = 6 pm). Measurements were taken on a sample with 
13x13 = 169 features, resulting in increments in the cavity 
separation parameter d of d step =118 nm. 

Once the two mirrors are brought close to each other 
within tens of microns, they are made parallel using a Thor- 
labs kinematic mount and the Fabry-Perot fringes present in 
the planar part of the mirrors with an accuracy of 200 prad. 
The cavities are placed close to one of the rotation axes of 
the mount, since this configuration enables us to tune the rel¬ 
ative angle between two adjacent coupled cavities without 
inducing significant cavity length changes (see Fig. SI b)). 
The setup is then initialised to a cavity length of L ~ 3 pm 
by obtaining the free spectral range with white light trans¬ 
mission (see Fig. SI c)). The mode properties of individual 
cavities are now accessible with a fixed wavelength laser 
transmission experiment. By sweeping the cavity length 
with a piezo-microactuator the position of modes is mea¬ 
sured relative to a defined ground state, thus enabling the 
conversion between cavity length and spectral position (see 
Fig. SI d)). The length scale for this sweep is obtained by 
traversing successive ground states, for which the cavity 
length differs by |. After spatially selecting a single feature, 
each mode is identified by counting its nodes in x and y 
direction. Another way of obtaining this spectral and spatial 
information is used in the strong coupling demonstration, 
where we used the fluorescence of a dense layer of quantum 
dots to populate the cavity modes, which is conceptually 
similar to positioning a broad-linewidth, red lamp between 
the mirrors (for more detailed experimental description see 
Supplement 2). 


In the plano-concave configuration a tunable, open-access 
microcavity consists of two freely positionable, opposing 
mirrors, one flat and the other concave. Starting from flat 
Silica substrates from UQG Optics, we pattern one of them 
with a Focused Ion Beam (FIB) [30] in order to produce con¬ 
cave templates. Our fabrication method produces the desired 


3. Mode formation theory and modelling 

The solutions of the paraxially approximated Helmholtz 
equation for cavities in the plano-concave geometry are the 
well-known Hermite-Gauss modes M q ^ n (x, y, z ). A state 
is defined by three natural numbers q , m and n , where q 
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Figure 2: Mode hybridisation following the formation rule ((m,^) + , ((2m, n), (2m + 1,«)) for merging fea¬ 

tures. 

a) Real part of wavefunction showing the hybridisation of the first bonding and anti-bonding states (upper two 
rows), and microscope images of the respective coupled cavity feature (lower row), b) Experimental results (dots) 
obtained by laser transmission experiment and theoretical values (dashed lines) derived with modemixing formal¬ 
ism. Theoretical mode positions are shifted towards bigger d values by 350nm (see text), c) Experimentally re¬ 
solved polarisation splitting of first bonding state, leaving two peaks with mutually orthogonal polarisation (inset). 
A P m in denotes linewidth resolution limit, A P max is calculated with formula (3) from [31] for £ = 1. d) Theoreti¬ 
cal finesse values for first bonding and anti-bonding states obtained with modemixing formalism. The dips corre¬ 
spond to resonances with higher transverse modes of same (m,ri) parity for lower q values (as reported in [37]). 


denotes the longitudinal and m, n the transverse excitation 
quantum numbers, counting the number of nodes of the 
electric field in the respective planes. For an azimuthally 
invariant cavity those modes are (m + n-h l)-fold degenerate 
but slight deviations lift the degeneracy and produce energy 
differences between the states. Each of these states has a 
further twofold polarisation degeneracy. 

In the modemixing formalism introduced by Kleckner 
et al. [38] solutions of the field distribution for a real cavity 
environment are found by solving the eigenvalue equation 
M l^i) = Yi\ x i*i), where M is the modemixing matrix which 
describes coupling between modes upon a single round trip 
within the cavity. The eigenvectors of M are the electric 
field amplitudes of stable cavity modes (T^), where the 
eigenvalue Ji corresponds to the scaling factor the ampli¬ 
tudes undergo for each round trip in the cavity and | V E ; ) is 
a linear superposition of initial basis states | \f/j). A natural 
choice for |i//}) is given by the normalised H q ^ n {x,y,z) 
forming an orthonormal basis in the transverse plane [39]. 


The spatial aspects of the cavity lead to the construction 
of M, which is obtained by modulating the pairwise over¬ 
lap of all basis states with the phase accumulated relative 
to a given plane due to the shape of the concave mirror 
A (x,y) (see Supplement 1). To introduce further aspects of 
the real cavity environment such as the angle between the 
mirrors, we deform A(x,y) by applying rotation matrices 
inducing a tilt by angle 0 around the v axis (Fig. SI b)). It 
can be shown that the paraxial approximation holds for a 
single cavity with the parameters that define the basis states 
(L ~ 3 pm and R = 6 pm). The same is true for the studied 
coupled cavity modes, since their waist is larger than in the 
single cavity case, thus guaranteeing a smaller divergence. 

Now the full spectrum is obtained by looping over the 
cavity separation d in small increments. At each d value 
the spectrum of eigenvalues and eigenvectors of the cor¬ 
responding modemixing matrix M is analysed, leading to 
the identification of the low energy modes in question. The 
spectral position (finesse) is obtained by analysing the phase 
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(magnitude) of the eigenvalue Yi and the spatial wavefunc- 
tion follows from the composition of the eigenvector | 
(Supplement 3 contains further information on this process). 


4. Results 

4.1. Mode hybridisation 

Figure 2 shows the mode hybridisation process that links the 
cavity states while merging. The modes shown correspond 
to a longitudinal quantum number q = 9 or a cavity length 
around L ~ 3 pm. The evolution of the modes from fully 
separated to merged is given by the rules (m,n) + (2m, n) 
and (m,n)_ —>> (2m + l,n). Here the subscripts + and — de¬ 
note the in and out-of-phase relation between the two bare 
modes. The rule is illustrated in Fig. 2 a), where the real part 
of the wavefunction for the first bonding and anti-bonding 
state is plotted alongside microscope images of the coupled 
cavity feature. Fig. 2 b) depicts experimental and theoretical 
spectral information of the first six modes relative to the 
ground state m + n = 0. The premise that fully separated 
( d = 10 pm) and merged (d = 0 pm) cavities display the 
same mode structure is found to be only approximately cor¬ 
rect for our system due to the long distance influence of 
the added shapes zq (x — d/2, y) and zq (x + d/2, y) on each 
other. For separations d oo ? the mode structure coincides. 
For d = 0 pm, the effective radius of curvature measured via 
the transverse mode splitting A*L is given by R = — - tT v 

sin 2 ( J 

which gives a value of 6 pm d= 0.1 pm in excellent agree¬ 
ment with the fabrication method input parameters. Starting 
with a cavity separation of d= 10 pm a pronounced splitting 
between the modes (1,0) and (0,1) is observed. The in- and 
out-of phase modes of the two branches stay degenerate 
down to d = 8.1 pm. For separation d below 8.1 pm the de¬ 
generacy is lifted and the coupling between bonding and 
anti-bonding state grows. For higher n states this process 
sets in later and the splitting is less pronounced. In between 
d = 6 pm and d = 8 pm, the mode energy decreases for 
higher m modes due to the decrease of the confining trans¬ 
verse potential. For smaller d the modes gradually evolve 
to the single deformed cavity case and display a gradually 
declining splitting with full recovery of the [m + n + l)-fold 
degeneracy for d = 0 pm. Experimental and theoretical val¬ 
ues show qualitative agreement with slight deviations in the 
horizontal position. This mismatch could be accounted for 
by shifting the theoretical curve in Fig. 2 b) and d) by 350 nm 
towards larger d values. The same shift has been applied to 
all theoretical graphs throughout this paper. We propose that 
this mismatch stems from the neglected penetration depth 
into the DBR (Ldbr ~ 900 nm) in the modemixing theory. 
Theoretically a decrease of the effective length of the cavity 
leads to a decreased mode waist on the curved mirror and to 
smaller inter-modal coupling, effectively shifting the curve 
towards smaller d values. With this compensation we ob¬ 
tain quantitatively well matched mode splittings across the 
whole d parameter range. 


In general each of these modes is twofold degenerate 
with respect to polarisation. Fig. 2 c) shows the lifting of 
this degeneracy for the first bonding state for intermediate 
d values between d = 7.0pm and d = 8.6pm, where the 
eccentricity of the concave shape is maximal, leaving two 
well distinguishable cross-polarised peaks (see inset of Fig. 
2c)). Below and above these values no polarisation induced 
splitting is observed within the linewidth resolution limit 
given by the coating limited finesse values of « 1000 (limit 
shown by A P m i n ). The evolution of the splitting with d tends 
to rule out a birefringence originating from the coating. In 
addition, we have checked over identical samples that the 
lower energy mode is always polarized perpendicularly to 
the coupling direction. We therefore attribute the splitting 
to nonparaxial field components along the longitudinal axis 
of the cavity as described in [31]. This effect is large in 
our structure due to the small radius of curvature compared 
to reference [31]. The maximal value of the splitting were 
taken from [31], and gives A P max by evaluating AP = 
for the eccentricity £ = 1. As seen in Fig. 2 c), this predic¬ 
tion is in good agreement with our measurement for highly 
elliptical cavities (d ~ 7.5 pm). At this position, the coupled 
cavity profile transitions from saddle shaped to elliptical. 
For both higher and lower d values, the splitting rapidly 
decreases below the linewidth resolution limit. 

Fig. 2 d) depicts finesse values of the first bonding 
and anti-bonding state. By application of the modemixing 
formalism each mode was assigned a finesse by equating 
Fi « y ^2 > where is the eigenvalue corresponding to the 

eigenmode in question [38]. The experimental limit is given 
by the dashed line and shows that the coating is the limiting 
factor for our system. The higher valued diffraction limits 
show a decrease for intermediary d values with several sharp 
dips in between before reaching maximal values for d above 
7.5 pm. To this end we capped the maximal finesse values at 
F = 10 5 corresponding to mirror losses of 30 ppm because 
of the finite numerical accuracy given by the modemixing 
approach. The general trend of the values shown is a result 
of the long distance influence of the two added concave 
shapes causing larger feature sizes on the high d end thus 
reducing diffraction losses. This effect vanishes for larger d 
values where the finesse returns to the values reported for 
d — 0 pm. The intermediary dips occur at positions where 
the respective mode is in resonance with a higher transverse 
mode of lower q which has the same parity. In this manner 
the bonding state (#, 0,0) has a dip at d = 3.85 pm where it 
is in resonance with the (#—1,2,2) state. This phenomenon 
has recently been reported for single cavities by Benedikter 
et al. [37]. Experimentally the resolution of these dips was 
not possible due to the relatively low coating-limited finesse 
(dashed line in Fig. 2 d)). 

Following the mode evolution for higher excited trans¬ 
verse states the modemixing formalism predicts avoided 
crossings between modes of same (m,n) parity. One of 
these anticrossings occurs between state (1,2) and (5,0) in 
single cavity notation around d ~ 6.31 pm with a splitting 
of h£l = 0.26meV. Fig. 3 shows the analysis of this phe¬ 
nomenon, depicting the spatial wavefunctions (a), the mode 
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Figure 3: Anti-crossing of higher transverse cavity modes 
with same (m,ri) parity, (1,2) and (5,0), obtained by 
modemixing analysis as the cavity separation is tuned 
across the resonance position, a) Spatial wavefunction show¬ 
ing gradual interchange of the mode identity, b) Spec¬ 
tral position (left ordinate) and finesse (right ordinate) 
demonstrating the avoided level crossing. The third red 
line corresponds to mode (3,1) in single cavity notation. 


positions (b, left ordinate) and the mode finesses (b, right 
ordinate) as obtained with the modemixing formalism by 
looping over d in 2.5 nm increments. The gradual exchange 
of the mode identity given by its finesse and spatial wave- 
function is accompanied by an avoided level crossing in the 
spectral domain. 


4.2. In-situ control of the coupling strength 

The wide length tunability of open-access microcavities 
allows a degree of in-situ control over the mode overlap 
and therefore inter-cavity coupling strength. In this set of 
experiments, the cavity length is scanned over the full sta¬ 
bility range to extract spectra for coupled modes as a func¬ 
tion of longitudinal mode number. Fig. 4 a) and b) present 
these logarithmic spectra stacked on top of each other for 
d = 7.77 pm and d = 7.41 pm. The bonding state mode posi¬ 
tion is used as a reference point set to 0 nm. Each spectrum 
is labelled with its respective cavity length for this ground 
state on the right side of Fig. 4 a) and b). The smallest op¬ 
tical cavity length achieved is 1.4 pm which is the smallest 
length attainable without putting the mirrors into contact 


Cavity 
length 

- 3.64 nm 
• 3.32 nm 
» 3 nm 

p 2.68 nm 
■ 2.36 nm 
1 2.04 nm 
P 1.72 nm 
b 1.4 nm 

0 50 100 ~ 

Relative cavity length (nm) 

Figure 4: Logarithmic transmission spectra for different cav¬ 
ity lengths (equivalent to longitudinal mode number q) for 
a feature separation of a) d = 7.77 pm and b) d = 7.41 pm. 
The experimental curves are shifted with respect to each 
other for clarity. The bonding (B) and anti-bonding (AB) 
states are highlighted with the black dashed lines. The abso¬ 
lute cavity length at the B state is given on the right hand 
side of the figure. 



(including the DBR penetration depth). Above L « 4 pm, 
the cavity mode linewidth increases significantly because of 
diffraction losses [38]. For L =3 pm the mode spacing corre¬ 
sponds to the values shown in Fig. 2 b) for the respective d 
parameter. In both Fig. 4 a) and b), the bonding (B) and anti¬ 
bonding (AB) states are highlighted with black dashed lines. 
The splitting between these two states is increasing with the 
absolute cavity length. This is consistent with the fact that 
the mode waists, and therefore the mode overlap, is increas¬ 
ing for larger length leading to an increased inter-modal 
coupling. A similar behaviour can be observed for higher 
excited states as well. For the cavity with d = 7.77 pm (re¬ 
spectively d = 7.41pm), the splitting is then tuned from 
2.7nm to 8.6nm (respectively 17.7nm to 28.3nm) by steps 
of 0.8 nm (respectively 1.7 nm). 


4.3. Strong coupling 

In order to prove that the cavities are strongly coupled and 
that the spectral splittings are not just an effect of the cav- 
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Figure 5: Mode splitting of ground (0,0) and first ex¬ 
cited (1,0) (bonding and anti-bonding) state for coupled 
cavity with d = 7.45 pm showing anticrossing behaviour 
when tuning the angle between the two mirrors (see SI 
b)). The dots correspond to experimental results with 
linewidths given by the error-bars and agree well with the 
theory (lines). The small insets depict the spatial wavefunc- 
tion for the two modes as observed experimentally (grey 
colourscale) and predicted theoretically via modemixing 
(blue-red colourscale). The large inset shows a selection of 
the original spectra; each mode is polarisation split due to 
the eccentricity of the cavity (comp, with Fig. 2 c)). 


ity shape, we demonstrate the characteristic anti-crossing 
behaviour between the ground state and first excited state. 
This experiment is performed simply by changing the angle 
between the two mirrors, thus varying the relative length 
of the two cavities. This experiment was performed for a 
separation parameter d = 7.45 pm and for the angular range 
—0.25° < 0 < 0.15° (see Fig. Sib) for visualisation). For 
each angle, the length of the cavity was adjusted to fix the 
ground state n + m = 0 at the spectral position A = 640 nm. 
Each mode was split by 0.3 nm into two orthogonally po¬ 
larised peaks due to the cavity eccentricity and was fitted 
with a double Lorentzian function (see the large inset in 
Fig. 5). The width and the position of the two peaks showed 
the same behaviour and contain the same information, which 
is plotted in Fig. 5 as a function of the angle 0. The dots 
denote the mean value of the fitted peaks and the error-bars 
give their respective linewidths. The continuous line is the 
result obtained with the modemixing formalism with rota¬ 
tion matrices M x (0) applied to modulate the mirror shape 
A(v,y). The anticrossing behaviour is clearly observed. For 
Fig. 5, the strength of the coupling, defined by the ratio of 
the mode splitting Q and the linewidth of the modes T, gives 
a value of / = p = 2, demonstrating the strong coupling 
character of the system. For large separation d, the strong 
coupling regime vanishes for f — 1 which corresponds, in 
our case to a separation of d = 8.1pm. Further pictorial 
confirmation of the accurate description by the modemixing 
formalism is given by the spatial wavefunctions obtained 
theoretically (red-blue colourscale insets) and experimen¬ 


tally (grey colourscale insets). The images show clearly the 
hybridisation between individual cavity modes for tuning 
through the avoided crossing. The same behaviour is seen 
for higher excited supermodes, which show less pronounced 
angular sensitivity the higher the excitation number m in the 
elongated direction is. 


5. Conclusions and Outlook 

In this work we have demonstrated control over cou¬ 
pled open-access microcavities with mode volumes around 
10 x A 3 . We presented experimental results from coupled 
cavities in which the inter-cavity spacing is varied and com¬ 
pared the results with predictions obtained by a modemixing 
formalism. Furthermore, the strongly coupled character of 
the bonding and anti-bonding states is demonstrated by the 
mode anti-crossing. By making use of the degrees of free¬ 
dom which the open-access cavity provides in-situ (length 
and angular tuning), the coupling between individual cavity 
modes of the system can be controlled. 

This work shows the ability to engineer the mode struc¬ 
ture of coupled open-access cavities as a fundamental build¬ 
ing block for photonic devices. Such photonic molecules 
have direct implications for the study of coupled quantum 
systems. For example they are attractive for the study of 
polariton blockade effects and the creation of ultrafast single 
photon sources. In addition the control of the photonic po¬ 
tential offers the prospect for engineering the Hamiltonian 
of coupled Bose-Einstein-Condensates [40]. Such studies 
could be applied to both exciton-polariton and photonic 
systems. 
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Supplemental Materials 

The shape of the concave feature forming one side of 
the coupled cavity is given analytically. Additionally the 
details of the experimental configuration and devices used 
to record the presented data in the main text is discussed. 
Finally the technical details of the modemixing formalism 
and the entailed computational methods are exposed. 


SI. Chip design 


We chose the shape of a single cavity to be the isophase 
surface of a Gaussian mode [39]. We keep this formulation 
for radial displacements that extend well beyond the paraxial 
approximation. This approach provides us an analytical 
form as well as a smooth shape which is ideal to prevent 
any crack formation in the DBRs. The phase of a Gaussian 
cavity mode at a position (x,y,z) is given by: 

ww)=M ^^-"(£) +fe (si> 


where R(z) is the radius of curvature (RoC) for v = y = 0 
and is given by R(z) = z l + (y ) 2 , Zr = (ttwq)/A is the 


Rayleigh length and k is the mode wavevector given by 
k = 2n/X. In order to confine the light, the cavity length L 
has to be smaller than the radius of curvature R on axis. 
In this case, the waist of a Gaussian mode is given by 

Wq = ^ \J wor k, we have used a RoC of 
6 pm for a targeted cavity length in the middle of the sta¬ 
bility range, ie. L = 3 pm. With this set of parameters, the 
waist is of the order of 800 nm, which leads to a single cavity 

mode volume of 44 x .To obtain the surface coordi¬ 
nate Z(j) as a function of (x,y) of the isophase surface, one 
has to solve the following equation: 



Figure S1: Gaussian mode isophase. The dashed lines corre¬ 
spond to the isophase for two independent cavities separated 
by 8 pm, the solid blue line to the sum of each independent 
isophase (eq. S4) taking the normalisation factor into ac¬ 
count. The black solid line gives the patterned shape which 
includes only the last 400 nm of the coupled cavity shape. 


the cavities presented in this article, these two cases are de¬ 
limited by a separation d = 7.6pm. We use the shape given 
by the function A (x,y) as an input of our fabrication method. 
The lateral extension of the feature is dictated by its depth 
that we choose to be 400 nm. For a single cavity, a depth 
of 400nm corresponds to a diameter around 5.5 pm, which 
is 7 times larger than the expected waist. Therefore, above 
99.99% of the mode is contained within the feature. The 
normalization term in equation S3 enables the comparison 
between each feature since the respective fundamental mode 
energies will be close to each other. 


ty( x ,y,z$) = 0(0,0,L) =kL — atan 



(S2) 


The coupled cavity shape A (x,y) is constructed by adding 
two surfaces defined by the function z </> (x,y) separated by a 
distance d (blue dashed lines on Fig. SI a)) which leads to 
the equations: 


A (x,y)=Z(x,y)dL- (S3) 

Z(x, y) =Z(f>(x- d/2, y)-hZ(f>(x-h d/2, y) (S4) 

where is a normalization term that keeps the max¬ 
imum length of the coupled cavities constant. As long as 
there is a saddle point in between the two cavities, Z max « L 
and the normalization factor is close to unity. However, 
when the saddle point disappears, the normalization factor 
keeps a constant length for the new structure created. For 


S2. Experimental setup 

Initially the absolute cavity length is set to 3 pm using a 
white light transmission experiment (Bentham WLS100). 
The collected light is then focused on an Andor combined 
spectrograph/CCD with a 300 grooves/mm grating in order 
to record the transmission spectrum. 

We use laser transmission to study the spectral and spa¬ 
tial mode structure of the cavities. A continuous wave TOP- 
TIC A laser DL100 with a fixed wavelength of A = 637 nm 
(linewidth < 1 MHz) is focused with an aspheric lens with 
(focal length = 5 cm) onto the chip. In this way the spot 
size is always bigger than the cavity contour allowing mode 
matching even to the higher transverse modes discussed 
in the main text. The chip is then imaged onto a Andor 
Newton CCD or an avalanche photodiode from Thorlabs 
(ref: APD120A2/M). One of the mirrors is mounted onto a 
ring piezo-actuator from Piezomechanik in order to sweep 
the cavity length. The length scale for the cavity sweep 
is obtained by traversing through successive ground states 
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with a well known length separation of j . By spatially se¬ 
lecting specific cavities, we can record the transmission 
spectrum as a function of cavity length. With this method, 
each maximum transmission peak is associated with a spa¬ 
tial wavefunction imaged onto a CCD allowing us to count 
the number of nodes and anti-nodes along the x and y di¬ 
rection, and therefore, properly identify each mode. Since 
the fundamental mode energy changes significantly as a 
function of d, we do not have an absolute measurement of 
the cavity length. Therefore, each mode observed is referred 
to with its energy difference from the ground (bound) state. 

A different method of obtaining single cavity spectra is 
used in the strong coupling demonstration, where the an¬ 
gle between the two mirrors is varied. Here we embed an 
ensemble of core-shell CdSe/ZnS quantum dots (Lumidot 
640 from Sigma Aldrich) into a thin spin-coated PMMA 
film on the planar mirror. After initialising the setup to 
a similar cavity length as above (L ~ 3 pm) with a white 
light transmission experiment, a PicoQuant LDH-D-C-470 
laser with A = 470 nm and continuous wave mode is fo¬ 
cused at a single cavity with moderate power densities of 
~ 1000 W/cm 2 to excite the quantum dots and cause emis¬ 
sion into the cavity modes around A = 640 nm. By adapting 
the cavity length such that the m + n = 0 mode is positioned 
at A = 640 nm the relative mode spacing to higher excited 
states is determined with lorentzian fits to the mode peaks 
in the cavity emission spectrum, which is recorded with a 
1800 grooves/mm grating. 


bounds of 10“ 12 and 10“ 10 respectively running on a paral¬ 
lel pool with 4x8 workers each running at 2.6 GHz. 

We infer the mode wavelength from the phase A(pj of the 
eigenvalue by noticing that A ? = —, where A and L are 

^ 4 nL 

the wavelength and cavity length defining the original basis 
states. The mode finesse is obtained by equating Fi « > 

which corresponds to diffraction losses only [38]. This pro¬ 
cedure leads to spectral and spatial determination of the 
mode properties. The full spectrum of the different cavity 
separation states is obtained by looping over d in 250 nm in¬ 
crements with a finer mesh of 62.5 nm in the coupled regime 
for 6.8 pm < d < 8 pm. 


S3. Modemixing formalism 

In the modemixing formalism stable cavity modes are found 
by solving the eigenvalue problem = yffi, where 
M G Mat(n, n),i G 1,2,..., n and is an eigenvector with 
its corresponding eigenvalue Yi- The modemixing matrix 
M is constructed according to [38] eq. (4) to (6). Each 
eigenmode | T*;) is a linear superposition of initial basis 
states | y/j). A natural choice for the 11//}) is given by the 
normalised Hermite-Gauss basis states (jc, y, z) as de¬ 
fined in [39]. Here j G 1,2,..., n maps each dimension of the 
original modemixing matrix to a combination of transverse 
excitation numbers m, n. For our purposes a set of the first 
500 basis states H q , m , n (x,y,z) ordered by increasing m + n 
value is chosen. Due to the elongation of the cavity profile 
in one dimension we truncate the basis set at n + m= 15 
in the perpendicular direction allowing modes of order up 
to n + m = 38 in the elongated direction. For feature sepa¬ 
rations above d = 8 pm we change the origin of the basis 
set to co-align with the left feature. Care was taken that 
both methods produced the same results for an overlapping 
region 8 pm < d < 8.5 pm. The convergence of the process 
is evident by the fast decaying contributions of higher order 
modes in the composition of the eigenstates of M, leading 
to an accurate representation of the actual field distribution. 
The exact calculation of the matrix elements is facilitated 
by the symmetry of the system, allowing the integration 
to be done in one quadrant only. The method is based on 
MATFABs quad2d routine with absolute and relative error 
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